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Abstract 

With the global color symmetry model (GCM) being extended to finite chemi- 
cal potential, the density dependence of the bag constant, the total energy and the 
radius of a nucleon, as well as the quark condensate in nuclear matter are investi- 
gated. A maximal nuclear matter density for the existence of the bag with three 
quarks confined within is obtained. The calculated results indicate that, before the 
maximal density is reached, the bag constant, the total energy of a nucleon and the 
quark condensate decrease gradually, and the radius of a nucleon increases, with the 
increasing of the nuclear matter density. Nevertheless no sudden change emerges. 
As the maximal nuclear matter density is reached, a phase transition from nucleons 
to quarks takes place and the chiral symmetry is restored. 
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1 Introduction 



It is well known that the nucleonic and mesonic degrees of freedom play essential roles 
in describing properties of nuclear matter and finite nuclei. Meanwhile, nucleons bound 
in nuclear medium alter their properties from those in free space. Although how much 
those properties change is of fundamental interest in nuclear physics, it is still not clear 
now. Then, in the relativistic mean field theory (RMF) based on the a — u model, which 
is currently regarded as a standard model of nuclear physics, the effective nucleon mass 
in nuclear medium is determined in self-consistent iteration on the fields of nucleons and 
mesons (see, for example, Ref.fl]]). However, this model is valid only at the hadronic level, 
at which a nucleon is described as a point-like particle. Because of the importance of the 
substructure of nucleon observed in deep inelastic scattering experiments and predicted 
by quantum chromodynamics (QCD), it is imperative to employ the quark and gluon 
degrees of freedom to describe nuclear phenomena. Unfortunately, it is now very diffi- 
cult to make quantitative predictions with QCD at low and intermediate energy region. 
Therefore various phenomenological models based on the QCD assumptions, such as the 
bag models^, ||, quark-meson coupling (QMC) model[||], and so on, have been developed. 
In bag models, the fundamental quantities to characterize a nucleon is the bag constant, 
the radius and the energy of the bag. In the QMC model, nucleons are regarded as non- 
overlapping bags interacting through the exchange of mesons. The size and mass of a 
nucleon are represented by the radius and energy of the bag, respectively, too. In order 
to take the modification of the volume energy on the mass of a nucleon into account, the 
bag constant is also exploited in the QMC model. Originally, the bag constant held fixed 
as its free-space value in bag models 0], then a density dependent bag constant is intro- 
duced by a phenomenological dependence on the medium density |5|, || 0, H, U- However, a 
sophisticated QCD foundation of the density dependence of the bag constant, the radius 
and the mass of a nucleon is still lacking. 

Recently, an effective field theory of QCD, namely the global color symmetry model 
(GCM), has been developed [[It], |TT| , |I2| , p~3f| , and shown to be quite successful in describing 
hadron properties in free space(i.e., at temperature T = 0, chemical potential fi = 0). 
With the Dyson- Schwinger equation approach of QCD being extended to finite temper- 
ature, the deconfinement and chiral symmetry restoration, the 7r and p meson proper- 
ties, and a part of the baryon properties at finite-T and finite-/i have been investigated 
14) , |15| ) pL6|, |17|) [1S| , pTOfl . With the nontopological-soliton ansatz||, some parameters related 



to the nucleon structure, such as the bag constant, had also been introduced in the GCM 



at T = and \i = in a natural way[lC]. It seems that a QCD foundation for the bag 
models is proposed at the level of the GCM model. With the global color symmetry model 
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being extended to finite chemical potential /x, the dependence of the nucleon properties, 
such as the bag constant, the total energy and the radius of a nucleon on the nuclear 
medium density will be studied in this paper. 

The paper is organized as follows. In Section 2 we describe the formalism of the GCM 
model at finite chemical potential /x and the model of the relation between the chemical 
potential of quarks and the baryon density in nuclear matter. In Section 3 we represent 
the calculation and the obtained results of the bag constant, the bag radius and the bag 
energy, as well as the quark condensate as functions of the nuclear matter density. It 
contains also discussions on the results. In Section 4, a brief summary and some remarks 
are given. 



2 Formalism 



It has been known that the action of the GCM in free space (i.e., at chemical potential 
/x = 0) for the zero mass quark is given [|T0|| in the Euclidean space as 



S 



<Px(Tyq(x) 7 • 85(x - y) + A°B°(x, y) q(y) + / ef xef y 



B d (x,y)B%y,x) 
2g*D(x-y) 



where B e (x,y) is a bilocal field and A e are matrices of the Fierz transformation among 
the spin, color, and flavor spaces of the quarks. Extending it to finite chemical potential 
/x, we have the action of the GCM at finite chemical potential as 



d A xd A yq{x) 7 • 85 (x — y) — /174 + A 9 B e (x,y) q(y) + I d A xd A y 



S = 

from which the generating functional is defined as 

Z[f], v } = J DqDqDB e e [ - s+7lq+ ^ ] 
After integrating the quark fields, we obtain 

S = -Trlog [7 ■ 85(x -y)- /x 74 + A"B° , I d^xd^y 



B%x,y)B 9 {y,x) 
2g*D{x - y) 



4 j4 B B (x,y)B e (y,x) 



2gW(x - y) 



Generally, the bilocal field B (x,y) can be written as [13 



B%x, y) = Bl{x, y) + £ ^^(£±2) 



(2) 



where B®(x, y) = Bq(x — y) is the vacuum configuration of the bilocal field. In the lowest 
order approximation with only the Goldstone boson being taken into account, the §\ 
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includes a and 7r mesons, and fi (i = a, ix) stands for the decay constant of a, it mesons, 
respectively. The vacuum configuration can be determined by the saddle-point condition 
-^jw = 0, and a equation of the translation invariant quark self-energy E(g, fi) is obtained 

as 



S(p,/x) 



— — g D(p-q)—y, 



With g M = (q, (g 4 + i/i)) being introduced, Eq. (3) can be rewritten as 



(3) 



d A q t a 
u9 2 D(x -y)—j v - 



:1» 



(2tt)^ v aJ 2 1 2 7 -g + S(g,/i) '"2 
Taking the conventional decomposition for the quark self-energy 

E(p, f j) = i[A(p)-l]'yp + VB(p), 



(4) 



where V = a + m • T75 with restriction a 2 + 7? 2 = 
solving the rainbow Dyson-Schwinger equations 



1, one can fix the self-energy S(p, \i) by 



~2 



[A(p)-l]jr 



B(p) 



8 /• dS 



3 
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(27T) 



■g 2 D(p-q) 



A{q)q-p 



A 2 (q)q 2 + B 2 (q) 
B(q) 



(5) 



(6) 



(2vr) 4 ^ ^ ^ J A 2 (q)q 2 + B 2 (q) ' 
Basing on the solution of the Dyson-Schwinger equations, one can determine the bilocal 
field, and fix further the GCM action. 

With a nontopological-soliton ansatz, the action of the bag models at finite chemical 
potential \l in terms of the GCM can be given by extending the formalism proposed in 
Ref.0 as 



s b = 5^9i[*7 -d + Wo- a(a(x) - m{x) ■ fj 5 ))qi + S(a,7c,fi) 



where S(a, tt, //) includes only a and 7r mesons and reads 



f 2 



J 7T , 



3{a,ir,fi) = / [4H<W + ^{d^f - V(a,7r)]d*z + 



with 
V(a,7T 



(7) 



127T 2 

(2tt) 4 Jo 



s'ds' I log 



A 2 (s')s' + (a 2 + 7T 2 )B 2 (s') 



(a 2 + 7T 2 



1)5 2 ( S ') 



/4 2 0'V + £ 2 0' 



A 2 (s')s' + 5 2 (s') 

Analyzing the stationary property of the bag and differentiating Eq. (7), one has equations 
for the quarks and mesons 



[iy-d + /ry - a(a(x) - m(x) ■ 775)] <?; = 



(9) 



4 



oa(x) 

dS B , . 

oir(x) 

The quark field and a, ir meson fields in symmetric nuclear matter can be determined 
by solving the Eqs. (9-11) self-consistently. As a consequence, the corresponding energies 
can be obtained. From the restriction on the quark self-energy (Eq. (4)), it is apparent 
that the meson fields corresponding to the vacuum configuration can be simply taken 
as a — 1, 7T = under mean-field approximation. In light of the nontopological-soliton 
ansatz||, |I0| , one can take the meson fields as a = and 7r = inside a bag (i.e., in 
a nucleon). For the quarks in a bag, Eq. (9) can thus be rewritten in the mean-field 
approximation as 

[z 7 • d + /z 7 o] qi (x) = (12) 

The lowest total energy of a single quark with respect to the radius R of the bag is given 

as 

*(R) = (13) 

where ujq = 2.04. This is consistent with the results obtained in Refs@, ||, |J. And the 
bag constant B is obtained as 



12vr 2 f°° . , . f 



A 2 (s')s' + B' 
A 2 (s')s' 



B 



2/ J\ 



A 2 (s')s' + B 2 (s>) 



(14) 



with s f = p 2 . With the correction from the motion of center-of-mass, the zero-point effect 
and the color-electronic and color-magnetic interactions being taken into account, the 
total energy of a bag is given as 

E = 3e(R) + t*#B - f = + ^B , (15) 

where Z /R denotes the corrections of the motion of center-of-mass, zero-point energy 
and other effects. 

Just as the same as that in Ref . (T^j , the bag is identified as a nucleon in the present 
work. It satisfies then the equilibrium condition 

dE(R) 
dR 

From this condition, we get 
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where a = 3uj — Z . As a consequence, Eq. (15) can be rewritten as 



3 V a ) 



(17) 



It is apparent that, with the solutions of Dyson-Schwinger equations (Eqs. (5) and 
(6)) being taken as the input for Eqs. (14), (16) and (17), the properties of nucleons (i.e., 
bags) in nuclear matter can be obtained. Because the quark condensate has commonly 
been taken as the order parameter to characterize the phase transition for chiral symmetry 



restoration |20], |2T|, we evaluate also the variation of the quark condensate against the 
density of nuclear matter 

B(s') 



< qq >=<: o(0)g(0) :>= TrS(x,x) 
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(2*0 



s'ds'— 

A 2 (s')s' + B 2 (s') 



In the practical calculation, since the knowledge about the exact behavior of g 2 and 
D(p — q) in low energy region is still lacking, one has to take some approximations or 
phenomenological form to solve the Dyson-Schwinger equations. For simplicity, we adopt 
the infrared dominative formf22|, |T0| 



g 2 D(p-q) = —rj 2 5(p~q) 
Id 



(19) 



where rj is a energy-scale parameter and can be fixed by experiment data of mesons. 
Although this form does not include the contribution from the ultraviolet energy region, 
it maintains the main property of QCD in the low energy region. With Eqs. (5), (6) and 
(19), one has 



A(p) 
A(p) 



1 



\ + 2 4 

v 



1/2- 



B(p) 
B(p) 



(ji — Ap 



0. 



for 



for 



~2 V 

p <T 

P > T 



(20a) 
(206) 



In order to investigate the dependence of nucleon properties on the nuclear matter 
density p explicitly, we must transfer the above obtained /^-dependence to that of the p- 
dependence. Because of the fermionic properties of quarks, the bags can be approximately 
regarded as a Fermi-Dirac systems [p3j| with the fermion number density 



n 
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k F d 3 k 



(2vr) : 



(21) 



where g is the degenerate factor, and kp is the Fermi momentum. For the quarks in a 
nucleon, g is 12, and 



6 



where m q is the current quark mass. In case of zero current quark mass, one has kp = p. 
In the lowest order approximation, we get the quark number density as 

n q = — 2 ^. (22) 

Considering the fact that the baryon number for a quark is 1/3, one can take the 
relation between the nuclear matter density and the chemical potential as 

'- = T = 3^ 3 ' < 23 > 

Combining Eqs. (14), (16-18), (20) and (23), we can obtain the dependence of the bag 
constant £>, the total energy E and the radius R of a nucleon (i.e., those of a bag) and 
the quark condensate on the nuclear matter density p B . 



3 Calculation and Results 

By calibrating the nucleon mass Mq = 939 MeV as the total energy of a bag and radius 
R = 0.8 fm in free space (i.e., p = 0, p = 0), we get the energy-scale 77 = 1.220 GeV, 
Z = 3.303. Such best fitted energy-scale 77 fits well the value 1.37 GeV, which was 
fixed by a good description of 7r and p meson masses [HI, |l7l , and is much more close to 



the Bjorken-scale 1.0 GeV (see Ref. [23] and the references therein). The obtained Zq is 



larger than the originally fitted value 1.84fl2"4|. However what we refer to here includes 



all the effects but not only the zero-point energy. Meanwhile other investigations (see for 
example Ref.|25|]) have shown that the zero-point energy parameter can be larger than 
1.84, even though the other effects are taken into account separately. With the above 
parameters 77, Z and Eqs. (17) and (18), we get at first the bag constant, and the quark 
condensate in free space as B = (172 MeV) 4 , and < qq > = —(132 MeV) 3 . It is evident 
that these obtained values B and < qq > are quite close to the results given in Ref. 0] 
and Ref. |2(|, respectively. 

By varying the chemical potential /x, we obtain the relation between the nuclear matter 
density and the chemical potential of the quarks, which exhibits the same monotonousness 
apparently. Furthermore, we get the variation behavior of the ratio of bag constant, the 
nucleon radius, the total energy of the bag and the quark condensate in nuclear matter 
to the corresponding value in free space against the nuclear matter density. The results 
are illustrated in Figs. 1-4, respectively. 

Looking over the figures, one may easily realize that, as the density of nuclear matter 
increases, the bag constant, the total energy of the bag and the quark condensate decrease 
monotonously. Meanwhile, the radius of a nucleon increases. When the nuclear matter 
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Figure 1: Calculated ratio between the bag constant in nuclear matter and that in free 
space as a function of the nuclear matter density 



20 
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Figure 2: Calculated ratio between the radius of nucleon in nuclear matter and that in 
free space as a function of the nuclear matter density 
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Figure 3: Calculated ratio between the total energy of a bag in nuclear matter and that 
in free space as a function of the nuclear matter density 




Figure 4: Calculated ratio between the local quark condensate in nuclear matter and that 
in free space as a function of the nuclear matter density 
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density reaches the value larger than 12 times the normal nuclear matter density (referred 
to as po), the bag constant, the total energy of a bag and the quark condensate vanish 
simultaneously, and the radius of a nucleon becomes infinite. Such behaviors indicate 
that nucleons can no longer exist as bags consisting of quarks. Since the density 12p 
is larger than the maximal average hadron density of neutron stars, which is commonly 
believed to be about 10po, the presently obtained results show that the nuclear matter 
changes to quark matter, i.e., the phase transition from hadrons to quarks happens, as 
the density of nuclear matter gets beyond the maximal average density of neutron stars. 
It provides a clue that there may be quark matter or hybrid matter with hadrons and 
quarks in the center part of neutron stars. This is something similar to the result given 
in Ref.p6[. On the other hand, the quark condensate < qq > has been regarded to be 
a manifestation to identify the chiral symmetry breaking. The gradual decrease of the 
quark condensate indicates that the chiral symmetry is restored gradually as the nuclear 
matter density increases. When the nuclear matter becomes quark matter, the chiral 
symmetry is restored completely since the quark condensate vanishes. It indicates that 
the quark decoupling phase transition and the chiral phase transition may happen at the 
same nuclear matter density. It is also worth mentioning that the increase of the nucleon 
radius induces naturally the swell of nucleons, which is believed to be essential to the 
EMC effect (see for example Ref.@). 

To show the function of the bag in nuclear matter and the process of the phase 
transition, we evaluate the critical radius of the bag in nuclear matter and the variation 
feature of the pressure of the bag against the nuclear matter density. In a classical point 
of view, the relation between the critical radius of the bag r c and the nuclear matter 
density p reads 



From the thermodynamics of quark system |27|], the pressure in a bag with quarks 



V 



/ ndp — B . 

Jo 



The numerical results of the critical radius and the pressure of the bag are displayed in 
Figs. 5 and 6, respectively. It is evident that, as the nuclear matter is dilute, the pressure 
of the bag is negative so that the bag exists definitely and three quarks are bound in 
a bag to form a nucleon. Furthermore, the bags exist separately from each other, since 
the radius of a bag is smaller than the critical radius r c . When the density of nuclear 
matter increases, the bags get close to each other. As the nuclear matter density is about 
2 times the normal density of nuclear matter, the nucleons begin to overlap with each 
other. However, every nucleon prefers to exist as an independent bag since the pressure 
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Figure 5: Calculated critical radius of a bag in nuclear matter with respect to the nuclear 
matter density 

of the bag is still negative. As the density p > 3p , the pressure changes to positive. It 
is definite that such a positive pressure enhances the overlapping among nucleons. Then 
quark matter appears as all the nucleons overlap with each other completely as the nuclear 
matter density p > 12p . As a consequence, the phase transition from hadrons to quarks 
takes place. On the other hand, taking the same way as that in determining the relation 
between nuclear matter density and quark chemical potential, and considering the fact 
that the degeneracy of nucleons is g — 4, one can easily get that p B = -^sP B - Comparing 
this result with Eq. (23), one can know that the chemical potential of quarks is just the 
same as that of nucleons. Such a equivalence means that the phase transition can happen. 

Figures 1 and 2 show also that, before the phase transition takes place, the changing 
characteristics of the bag constant and the radius of a nucleon coincide with the results 
obtained in the QMC model with a phenomenological dependence on the scalar meson 
field being involved^, ||, |7|, ||. Such a behavior indicates that, dealing the bag constant and 
the bag radius with a phenomenological dependence on the medium density in the QMC 
model is reasonable. Nevertheless, Fig. 3 shows that the change of the total energy of a 
bag is much less rapid than that of nucleon mass obtained in QMC and other frameworks. 
Tracking down the source of the discrepancy, one knows that the total energy of a bag is 
not identical to the mass of a nucleon, but with a relation M = \JE 2 — < p > 2 . Taking 
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Figure 6: Calculated pressure of the bag with respect to the nuclear matter density 

< p > 2 = [i 2 B approximately, we get that, at the normal density of nuclear matter, the 
mass of a nucleon decreases to about 85% of the value in free space. This is consistent 
with the result determined self-consistently in the derivative scaler coupling model of 
the RMF||. It means then the effective mass of nucleons in RMF calculations is also 
reasonable. Such agreements show that present model can reproduce the quantities and 
their variation characteristics obtained in the QMC model and the RMF qualitatively. 

4 Summary and Remarks 

In summary we have investigated the density dependence of the bag constant of nucleons, 
the nucleon radius and the total energy of the bag as well as the quark condensate in 
nuclear matter in the global color symmetry model, an effective field theory model of 
QCD. A maximal density of nuclear matter, which is larger than 12 times the normal 
nucleon density, for the existence of the bag of quarks is obtained. The calculated results 
indicate that the bag constant, the total energy of the bag and the quark condensate 
decreases with the increasing of the nuclear matter density before the maximal density is 
reached. Meanwhile the size of nucleons swells. As the maximal density is reached, a phase 
transition from nucleons to quark-gluons takes place and the chiral symmetry is restored. 
Furthermore, the presently obtained changing features agree with those obtained in QMC 
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and RMF quite well. In this sense, it provides a clue of the QCD foundation to the QMC 
and RMF with a simple effective field model of QCD. 

In the present calculation, the g 2 D{p — q) is taken to be proportional to a 5— function. 
However, the detailed effects of the running coupling constant, the gluon propagator 
D(p — q) and the other degrees of freedom on the changing feature have not yet been 
included. Especially, since the meson fields were taken to be a = and ix = in the 
bag with respect to those of the vacuum configuration a — 1, 7r = 0, the self-consistent 
interaction and adjustment between the quarks and the meson fields have not yet been 
taken into account. Meanwhile, the relation between the chemical potential and the 
nuclear matter density was handled with a simple correspondence in statistical mechanics. 
It means that the present calculation is an approximation of the scheme. The obtained 
results are thus the preliminary ones. In principle, the quark energy (and field), the meson 
fields and the nuclear matter density can be determined by solving the set of differential- 
integral equations established from the action(Eq.(l)) of the system consistently. Such a 
sophisticated investigation is under progress. 
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